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NONLOCAL ELLIPTIC EQUATIONS IN BOUNDED DOMAINS: 

A SURVEY 

XAVIER ROS-OTON 


Abstract. In this paper we survey some results on the Dirichlet problem 

r Lu = / in n 
1 u = g in ]R"\r2 
for nonlocal operators of the form 

iit(x) = PV / {u{x) — u{x + y)^K{y)dy. 

We start from the very basics, proving existence of solutions, maximum principles, 
and constructing some useful barriers. Then, we focus on the regularity properties 
of solutions, both in the interior and on the boundary of the domain. 

In order to include some natural operators L in the regularity theory, we do 
not assume any regularity on the kernels. This leads to some interesting features 
that are purely nonlocal, in the sense that have no analogue for local equations. 

We hope that this survey will be useful for both novel and more experienced 
researchers in the field. 
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1. Introduction 

The aim of this paper is to survey some results on Dirichlet problems of the form 


Lu = f in Q 
u = g in 


( 1 . 1 ) 


Here, is any bounded domain in M”, and L is an elliptic integro-differential oper¬ 
ator of the form 


Lu{x) = FV / {u{x) — u{x + y)'^K{y)dy. 
The function K{y) > 0 is the kernel of the operatoiQH, and 


( 1 . 2 ) 


we assume 


K{y) = K{-y) and 


mm 


{\y\‘^A}K{y)dy < oo. 


Integro-differential problems of the form fll.ip arise naturally in the study of sto¬ 
chastic processes with jumps, and have been widely studied both in Probability and 
in Analysis and PDEs. We refer to the reader to the following works on: 

• Existence of solutions: Felsinger-Kassmann-Voigt ESI, Hoh-Jacob EH. and 
Barles-Imbert [3]; 

• Interior regularity of solutions: Bass-Levin [6], Kassmann [33], Caffarelli- 
Silvestre [niiiiiig, Barles-Chasseigne-Imbert [5], Kassmann-Mimica 
[35] . Schwab-Silvestre [5T], and Serra [55]; 

• Boundary regularity of solutions: Bogdan [8], the author and Serra 
l47] . Grubb [29l[30], Chen-Song [16], Barles-Chasseigne-Imbert [4], Bogdan- 
Grzywny-Ryznar [9], and Bogdan-Kumagai-Kwasnicki [lOj : 

• Other qualitative properties of solutions: Birkner-Lopez-Wakolbinger [7], 
Dipierro-Savin-Valdinoci EDI. Abatangelo [1], and Kulczycki 


In applications, this type of nonlocal problems appear in models involving “anoma¬ 
lous diffusions” (in which the underlying stochastic process is not given by Brownian 
motion), or in the presence of long-range interactions or forces. In particular, they 
appear in Physics [22l[58l|l2l[27l[23l[38] , Finance [53l[39ll2], Fluid dynamics [T^fTT] . 
Ecology [m [321 m], or Image processing 


The aim of this paper is to survey some results on existence, regularity, and 
qualitative properties of solutions to fll.ll) . More precisely, we will start from the 
very basics (showing existence and boundedness of solutions), to then focus on the 
regularity of solutions, both in the interior and on the boundary of the domain G. 


^The typical example is K{y) = Cn,s\y\~^~^‘^, which corresponds to L = (—A)®, the fractional 
Laplacian. As s f 1, it converges to the Laplacian — A. 

^For equations with x-dependence, the kernel is a function of two variables, K{x,y). Here, we 
assume that the operator L is translation invariant, and thus that the kernel K does not depend 


on X. 
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In order to include some natural operators L in the regularity theory, we do not 
assume any regularity on the kernel K{y). As we will see, there is an interesting 
relation between the regularity properties of solutions and the regularity of the 
kernels K{y). This is a purely nonlocal issue, in the sense that has no analogue in 
second order equations, as explained next. 

For linear elliptic second order equations of the form 

f = / infi 

1 u = g on do,, ^ ' 


the interior regularity properties of u depend only on the regularity of /. This is 
because, after an affine change of variables, this equation is just —Am = /. 

The nonlocal analogue of fll.311 is fll.ip . For this problem, the interior regularity 
of solutions depends on the regularity of / —as in fll.Sp — , but it also depends on 
the regularity of K[y) in the y-variable. Furthermore, if the kernel K is not regular, 
then the interior regularity of u will in addition depend on the regularity of ( 7 , on 
the boundary regularity of m, and even on the shape of fi, as we will see later on. 

Hence, for nonlocal equations, the class of linear and translation invariant oper¬ 
ators is much richer, and already presents several interesting features. The same 
type of issues appear when considering nonlinear equations, or operators with x- 
dependence. However, for the clarity of presentation, we will consider only linear 
and translation invariant equations. 

In most of the regularity results we will focus on two classes of kernels: we will 
assume s G ( 0 , 1 ) and either that 

\y\n+2s - - |^|n+ 2 s’ 0 < A < A, (1.4) 

or that 

= aeT(S'-), a>0. (1,5) 

In fact, in order to include operators like L = - 1 - • • • 

will allow a in fll.511 to be any nonnegative measure on 5'"'“^. The only necessary 
condition on the measure a is that it is not supported in a hyperplane —so that 
the operator L is not degenerate. A quantitative way to state this nondegeneracy 
condition is that, for some positive constants A and A, one has fgn-i da < A and 
fsn-i W ■ d\‘^^da{6) > A for all u G 5'"'“^; see [47] . 

The paper is organized as follows. In Section [2] we briefly explain the probabilistic 
interpretation of (II.ip . In Section |3] we show existence of weak solutions. In Section |4] 
we prove the maximum principle for such solutions. In Section [5] we construct 
some useful barriers and give an L°° estimate. In Section | 6 ] we discuss the interior 
regularity properties for the classes of kernels fll.4p and fll.Sp . In Section [7] we see 
what is the boundary regularity of solutions. Finally, in Section [ 8 ] we come back to 
the interior regularity of solutions to fll.ip . 
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2. Motivation and some preliminaries 


2.1. Levy processes. Integro-differential equations of the form (11. Ih arise naturally 
in the study of stochastic processes with jumps, and more precisely in Levy processes. 
A Levy process is a stochastic process with independent and stationary increments. 
Informally speaking, it represents the random motion of a particle whose successive 
displacements are independent and statistically identical over different time intervals 
of the same length. 

These processes extend the concept of Brownian motion, and were introduced a 
few years after Wiener gave the precise dehnition of the Brownian motion [ID]. Es¬ 
sentially, Levy processes are obtained when one relaxes the assumption of continuity 
of paths (which gives the Brownian motion) by the weaker assumption of stochastic 
continuity. 

By the Levy-Khintchine Formula, the inhnitesimal generator of any Levy processes 
is an operator of the formic 


= + / {n{x + y) - u{x) - y ■Vu{x)xBAy)]di^{y), 

ij j dRr. 


where u is the Levy measure, and satishes /jj^min{l, \y\‘^'^di'{y) < oo. When the 
process has no diffusion or drift part, this operator takes the form 


-Lu{x) = / {u{x + y)-u{x)-y■Vu{x)xBl{y)}du{y)■ 

JM.^ 

Furthermore, if one assumes the process to be symmetric, and the Levy measure to 
be absolutely continuous, then L can be written as fll.2p or, equivalently, 

Lu{x) = \ f {2u{x) - u{x + y)- u{x - y)]K{y)dy, 

with K[y) = K{—y). We will use this last expression for L throughout the paper. 

As an example, let C M” be any bounded domain, and let us consider a Levy 
process Xt, t > 0, starting at x G ff. Let u{x) be the expected hrst exit time, i.e., 
the expected time E[r], where r = inf{t >0 : W ^ 12} is the hrst time at which 
the particle exits the domain. Then, u{x) solves 

J Lu = 1 in 12 
} M = 0 in M”\12, 

where —L is the inhnitesimal generator of Xt. 

Recall that, when Xt is a Brownian motion, then L is the Laplace operator —A. 
In the context of integro-diherential equations. Levy processes plays the same role 
that Brownian motion play in the theory of second order equations. 

Another simple example is given by the following. Let us also consider a bounded 
domain 12, a process Xt starting at a; G 12, and the hrst time r at which the particle 


^We denote the infinitesimal generator — L in order to be consistent with our notation in (II.2|) . 
We always take L to be positive definite, so it is the analogue of —A. 
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exits the domain. Assume now that we have a payoff function g : M"" \ ^2 —)■ M, so 
that when the process Xt exits we get a payoff g{XT-). Then, the expected payoff 
u{x) := E[ 5 f(XT-)] solves the problem 

J Lu = 0 in n 

I u = g in 

The Dirichlet problem fll.ip arise when considering at the same time a running 
cost / and a final payoff g. 

2.2. Kernels with compact support. It is important to remark that when the 
kernel K{y) has compact support in a ball Bs (for some 5 > 0), then the Dirichlet 
problem is 

f Lu = f in fl 

\ u = g in 

This means that the exterior condition g has to be posed only in a neighborhood of 
dQ and not in the whole M”\r2. This turns out to be convenient in some applications; 
see for example | 22 l |56l |58l [TT] . 

However, from the analytical point of view, it is only the singularity of the kernel 
K(y) at the origin that plays a role in the regularity properties of solutions; see for 
example Section 14 in [T3] . 

2.3. Stable processes. A special class of Levy processes are the so-called stable 
processes. These are the processes which satisfy self-similarity properties, and they 
are also the ones appearing in the Generalized Central Limit Theorem; see |50] . 

The inhnitesimal generators of these processes are given by fll.2p - fll.5p . Note that 
the structural condition fll.Sp on the kernel K is equivalent to saying that the Levy 
measure is homogeneous. This is also equivalent to the fact that the operator L is 
scale invariant. 

A very natural stable process is the radially symmetric one. This means that 
K[y) = c\y\~'^~‘^^ and hence, up to a multiplicative constant, 

L = (-Ar. (2.1) 

Another very natural stable process is the one obtained by taking independent 
stable processes in each coordinate. That is, we consider = (xj^\xj"'^), 
where xj:^^ are 1-dimensional i.i.d. symmetric stable processes. The generator of 
this process X^ will be 

^ *1 -^ ( 2 - 2 ) 

and it corresponds to fll.Sp with the measure a being 2n delta functions on 
For example, when n = 2, one has a = (5(o,i) -l- <5(o,-i) + <5(i,o) + <5(_i,o)- 

While the operators fl2.ip and fl2.2p may look similar (they are the same when 
s = 1 ), they have quite different regularity properties, as we will see in the next 
sections. For example, while solutions to {—A)^u = 0 in D are C°° inside the 
domain, this may not be the case for fl 2 . 2 p . 
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Note also that the Fourier symbols of (I2.1jl and (I2.2p are and - 

respectively. Again, the first one is C°° outside the origin, while the second one is 
just In general, the Fourier symbol of any stable process of the form fll.Sp is 

^(0 = c f |^ 9\^^a{e)de. (2.3) 

The symbol A(^) is in general only outside the origin, but it would be C°° 
outside the origin whenever the function a is C°° on 5'"'“^. 

3. Existence oe solutions 

In this Section we explain briefly how to prove existence of weak solutions to fll.ip 
in a simple case. We hope that this will be useful for students and/or for readers 
which are not familiar with nonlocal operators. The more experienced reader should 
go to Felsinger-Kassmann-Voigt [25], where this is done in a much more general 
setting. 

As we will see, the existence (and uniqueness) of weak solutions to fll.ip - fll.2p 
follows from the Riesz representation theorem once one has the appropriate ingre¬ 
dients. 

The energy functional associated to the problem fll.ip is 

E{u) = - f f {u{x) — u{z)YK{z — x)dx dz — f fu. (3.1) 

4 J Jn 

The minimizer of £ among all functions with u = g m \ hi will be the unique 
weak solution of fll.ip . 

Notice that £{u) is dehned for all regular enough functions u which are bounded 
at infinity. When <7 = 0, then making the change of variables y = x + z the functional 
can be written as 

^(«) = T [ [ {u{x)-u{x + y)YK{y)dxdy- [ fu. (3.2) 

4 jR" Jr" Jn 

When g is not zero, the term f l5'(^) ~ 9(z)l^d^(z — x)dxdz could be infinite, 

and this is why in general one has to take fl3.ip . 

For simplicity, we will show existence of solutions for the case g = 0, and hence 
we can think on the energy functional fl3.2p . 

Let be the space of functions u G L^(M") satisfying 

['>Ahk-=\[ I {u{x)-u{x + y)YK{y)dxdy < oo, (3.3) 

JR" JR" 

and let 

X = {ue i/i^(M”) : u = 0 in \ fl} . 
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Essentially, the only assumption which is needed in order to prove existence of 
solutions is the Poincare inequality 


<C 


I 12 

\u{x) — u{x + y)\ K{y)dxdy 


(3.4) 


Jn /TOn 


for functions m = 0 in M” \ 

When the kernel K satishes fll.dp . then the seminorm in fl3.3l) is equivalent to 

|2 


[u\^s 


|■u(x) — u{x + y)\‘ 


dx dy. 


(3.5) 




^y^n-\-2s 

In this case, the Poincare inequality fl3.4p follows easily from the fractional Sobolev 
inequaht}0 in M” and Holder’s inequality in 

Once one has the Poincare inequality (13.dp . then it follows that the space X is a 
Hilbert space with the scalar product 


{v,w)k := - 


/ron 


{v{x) — v{x + y)) (w{x) — w{x + y)')K{y)dx dy. 

Then, the weak formulation of (II.ip is just 

{u, (p)k = / for all G X, (3.6) 

Jn 

and the existence and uniqueness of weak solution follows immediately from the 
Riesz representation theorem. 

For more general classes of kernels K{y), one has to show the Poincare inequality 
fl3.4l) (see Lemma 2.7 in [25]), and then the existence of solutions follows by the 
same argument above. 

Finally, notice that when u is regular enough then for all G X we have 
{u,ip) = - [ I {u{x) — u{z)){i.p{x) — i.p{z))K{x — z)dxdz 


iri ./ron 


= -PV 


Ml 


./ron 


1 

2 


Lu{x) (p{x)dx + 


x) — u{z)'j(p{x)K{x — z)dxdz+ 

u{z) — u{x)'j(p{z)K{x — z)dxdz 

Lu{z)ip{z)dz 


^PV 


$ri /ron 


= / LuLf, 

Jn 

where we used that K{y) = K{—y) and that (p = 0 in M"" \ 


^Different proofs of the fractional Sobolev inequality can be found in [^, [75], and [T5]. The 
proof given in [43] . which is due to Brezis, is really nice and simple. 
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This means that, when u is regular enough, the weak formulation (I3.6p reads as 



and thus it is equivalent to Lu = / in hi. 


for all G X, 


Remark 3.1. We showed in this Section how to prove the existence of weak solutions 
to fll.ip by variational methods. On the other hand, an alternative approach is to 
use Perron’s method to establish the existence of viscosity solutions; see [2US]. As 
we will see, when the right hand side / is Holder continuous then all solutions are 
classical solutions (in the sense that the operator L can be evaluated pointwise), and 
thus the notions of weak and viscosity solutions coincide; see [52] for more details. 


4. Comparison principle 

In this section we show the maximum principle and the comparison principle for 
weak solutions to fll.ip . 

As in the classical case of the Laplacian —A, the maximum principle essentially 
relies on the fact that Lu{xq) > 0 whenever u has a maximum at Xq. In case of 
local equations, this is true for any local maximum, while in nonlocal equations 
this is only true when u has a global maximum at xq- In this case, the inequality 
Lu{xo) > 0 follows simply from the expression 

Lu{xo) = \ j (2m(xo) - u{xo -\-y)- u{xo - y))K{y)dy 

^ Jw' 

and the fact that m(xo) > u{xq ± y) for any y G M” (which holds whenever u has a 
global maximum at xq). 

Notice also that when K{y) > 0 in M"", this argument already yields a strong 
maximum principle, since one has Lu{xo) > 0 unless u = ctt. 

The previous considerations work when u is regular enough, so that Lu can be 
evaluated pointwise. In case of weak solutions u to fll.ip . the proof of the maximum 
principle goes as follows. For simplicity, we do it in the case K{y) > 0, but the same 
could be done for more general kernels K{y) > 0. 

Proposition 4.1. Let L be any operator of the form fll.2l) . with K{y) >0 in M”. 
Let u he any weak solution to fll.ip . with f > 0 in Q and g > 0 in ML \Q. Then, 
u > 0 in Q. 

Proof. Recall that u is a weak solution of fll.ip if 

/ / {u{x) — u{z)'){ip{x) — y){z))K{z — x)dx dz = / fy) (4.1) 

for all y) G HKiM'^) with (^ = 0 in M” \ H. 

Write u = u'^ — u~ in H, where u'^ = max{M, 0}xn and u~ = max{—-u, 0}xn. We 
will take y> = u~, assume that u~ is not identically zero, and argue by contradiction. 
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Indeed, since / > 0 and > 0, then we clearly have that 

[ f^>0- 

Jn 

On the other hand, we have that 

(u{x) — u^z)) [(p{x) — ip{z)'jK{z — x)dxdz = 

{u{x) — u{z)) {u~{x) — u~{z))K[z — x)dx dz + 

+ 2 / dx {u{x) — g{z)')u~{x)K{z — x)dz. 

Jn Jn^ 

Moreover, (u'^{x) — W^^z)) (u~{x) — u~{z)'j < 0, and thus 

(^u{x) — u^z)) (u~{x) — u~{z)'jK{z — x)dx dz 


<— / (u {x) — u {z)yK(z — x)dxdz<0. 

Jn Jn 


(4.2) 



n Jn 



n J n 


Also, since g > 0 then 


Jn Jn<^ 
Therefore, we have shown that 


dx / {^u{x) — g{z))u {x)K(z — x)dz < 0. 



and this contradicts (I4.1l) - (l4.2p . 


(u{x) — u^z)) [(p{x) — ip{z)'jK{z — x)dx dz < 0, 


□ 


Of course, once we have the maximum principle, the comparison principle follows 
immediately. 

Corollary 4.2. Let L he any operator of the form (II.2p . with K{y) >0 m M”. Let 
ui and U 2 he weak solutions to 

j Lui = fI in kl , J Lu 2 = f 2 inVL 

\ ui = gi in M"'\0 \ U 2 = g 2 in M"'\0. 

Assume that fi > f 2 and gi> g 2 - Then, Ui > U 2 - 

Proof. Just apply Proposition 14.11 to u = Ui — U 2 . □ 


This allows us to use barriers, which in turn can be used to show that solutions 
u to fll.ip belong to L°°{fl) whenever / and g are bounded. This is what we do in 
the next Section. 
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5. Barriers and L°° bounds 
We provide in this Section an L°° estimate of the form 

||'w||l°°( 0) ^ ||fi'||L°°(R"\0) + C'||/||l°°(0) (5-1) 

for solutions to fll.ip . 

To do it, we assume hrst K{y) > {] in M"’, as in the previous Section. In this case, 
the construction of a barrier is quite simple. 

Lemma 5.1. Let L he an operator of the form fll.2p . with K{y) >0 in M”. Then, 
there exists a function w G (M”) such that 


Lw 

> 

1 

in 

n 

w 

> 

0 

in 


w 

< 

C 

in 

n. 


The constant C depends only on the kernel K and diam{Q). 

Once we have this barrier, the L°° bound (15.1 p follows from the comparison prin¬ 
ciple, as shown next. 

Corollary 5.2. Let L be any operator of the form (11. 2p . with K{y) > 0 in M"’. Let 
u be any weak solution of (II. ip . Then, 

where C is the constant in Lemma \5.1\ 

Proof. Let v{x) = || 5 f|| LjOO + 11/11 LjOO w{x), where w is given by Lemma [5.11 Then, we 
clearly have Lu < Lv in Q, and < n in M"" \ O. 

Thus, by the comparison principle, we have n < n in O. In particular, u < 
+ C||/||loo in O. 

Applying the same argument to (—n), we hnd that —u < ||( 7 || LjOO + C\\f\\ LjOO ^ and 
hence the result follows. □ 

We now construct the barrier. 

Proof of Lemma \5.1\ Let Br be any large enough ball such that O CC Br, and let 
T] G Cf°{BR) be such that 

0 < < 1 in M”, rj = 1 in O. (5.2) 

Then, for each a: G O we have ri{x) = maxRn rj, and thus 

2rj{x) — rj{x + y) — rj{x — y) > r]{x) — rj^x + y) >0. 

Hence, writing z = x + y, we have 
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where we have used that ri{x) — ri{z) = 1 in M” \ Bji. Now, we notice that 


^\Br 


K{x — z)dz = 


'^\{x+Br) 


K{y)dy > 


'^\B2R 


K{y)dy = c > 0 


for some positive constant c. 

Hence, we have Lr; > c > 0 in Taking w = -? 7 , we will have that 


Lw 

> 

1 

in 

n 

w 

> 

0 

in 


w 

< 

1 

c 

in 

n, 


as desired. 


□ 


Note that the previous proof works not only for all kernels satisfying K{y) >0 in 
”, but also for all kernels K satisfying 


^\Br 


K{y)dy > 0 


(5.3) 


for every ball In particular, the above construction works for all stable opera¬ 
tors fii.sp . 

Hence, we see that the desired barrier is quite easy to construct, and does not 
even depend much on the class of kernels we are considering. Essentially, thanks to 
the nonlocal character of the operator, any function rj as in fl5.2p works. 


Remark 5.3. For operators L that do not satisfy (I5.3p . one can still manage to prove 
an analogous result. Indeed, if (I5.3p does not hold, then K = Q a.e. outside a ball Br. 
This means that K have compact support. In this case, one can take any function 
r; > 0 which is strictly concave in a large ball Bm and zero outside Bm- Then, if M 
is large enough (so that f2-|-supp(it') C Bm and hence 2ri{x) — ri{x + y) — ri{x — y) > 0 
therein), one will have Ly > c > 0 in Q. 

Finally, to end this Section, we give an important explicit solution for the class 
of stable operators fll.5p . 

Lemma 5.4. Let L be any stable operator of the form fll.2p - fll.5p . Then, the func- 
tio"^ 


solves 


Uo{x) := (l - 

Luq = c in Bi 
Uq = 0 in MT\Bi 


for some positive constant c > 0. 


®Here z+ denotes the positive part of the number z, i.e., z+ = max{2,0}. 
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This explicit solution will be used in Section [7] to construct a subsolution and 
establish a Hopf Lemma for this class of operators. 

Lemma 15.41 was first established for (—A)® in dimension n = 1 by Getoor [2B] . 
and requires quite fine computations; see also the work of Dyda [21] . 

Once this is established in dimension n = 1, the result in dimension n and for 
general stable operators L follows by just writing 


Luq{x) = 7 / a{e) 
4 


2uo{x) — uo{x + t 9) — uo{x 

I 'j -1 lH-2s 



de, 


and using the result in dimension 1 for each direction 6. It is important to notice 
that the ID functions r Uo{x + t9) are exactly a rescaled version of (1 — |rp)^, 
which solves the equation in dimension n = 1. Using this, one gets that Luq = c in 
i?i, with c = Cg Jgn-i a{9)d9, and Cg depending only on s. 


6. Interior regularity 

For second order equations, the classical Schauder estimate for —An = / in i?i 
establishes that 

||w||c2+“(Bi/2) < C (||/||c“(Bi) + l|w|U°°(Si)) (6.1) 

whenever a > 0 is not an integer. Thus, it immediately follows from this estimate 
that solutions to the Dirichlet problem 

f —An = / in D 
u = g on d^l 

are inside D whenever / G G" and g is bounded. 

For nonlocal equations of order 2s, one would expect a similar estimate, in which 
the norm \\u\\c‘^s+a(^Bi/ 2 ) controlled by ||/||c“(Bi) for solutions to Lu = f in Bi. 
It turns out that, due to the nonlocality of the equation, the norm ||n||L°°(_Bi) in 
fl6.ll) has to be replaced by a global norm of n, i.e., a norm that controls n in the 
whole M". 

6.1. Regular kernels. For the fractional Laplacian L = (—A)^ —which corre¬ 
sponds to K{y) = in fll.2p — , this estimate reads as 

||m||c2»+“(Bi/2) a C (ll/llc^fBi) + ||m||l°°(R'‘)) , (6.2) 

and holds whenever a -|- 2s is not an integer; see for example [371 SS] • 

The same estimate fl6.2p holds for more general operators of the form fll.2p - fll.4p 
or fll.2p - fll.5D under the extra assumption that the kernels K{y) are G" outside the 
origin; see [SSI Corollary 1.2] and [HI Corollary 3.5] for more details. 

This means that, in this case of regular kernels, solutions to fll.ip are 
inside D whenever / G G“ and g is bounded. 
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However, for general operators (ll.2jl - (ll.4jl or (ll.2jl - (ll.5p (with no further regularity 
assumption on the kernel K), the estimate fl6.2p is not true anymore, and one needs 
a stronger norm of u in the right hand side, as explained next. 

6.2. Non-regular kernels. When the kernels K{y) are not regular outside the 
origin, one has the following estimate. 

Theorem 6.1 ([551117]). Let L be any operator of the form fll.2p with kernel K of 
the form either fll.ip or fll.5p . Let a > 0 he such that a + 2s is not an integer, and 
let u G L°°(ML) be any weak solution to Lu = f in Bi. Then, 

\\u\\c^s+a(^Bi/ 2 ) — ^ (ll/llc'“(Bi) + ||'*l||c“(R")) ) (6-3) 

for some constant C that depends only on n, s, and X, and A. 

It is important to remark that the previous estimate is valid also in case a = 0 
(in which the C" norm has to be replaced by the L°°)-, see Theorem 1.1 in [17| for 
more details. 

The proof of this estimate uses a rehned version of the blow-up and compactness 
argument hrst introduced by Serra in [51|. The estimate 06.31) was established in 
[55] for the class of kernels Ol.ip R. and in [IT] for the class 01.51) . 

With no further regularity assumption on the kernels K, the estimate 06.31) is 
sharp, in the sense that the norm ||u|lc“(R") can not be replaced by a weaker one. 
More precisely, one has the following. 

Proposition 6.2 ([551 ST]). Let s G (0,1), a G (0,s]. Then, for any small e > 0 
there exists an operator of the form 01. 21) - 01. 5p . and a solution u to Lu = 0 in Bi 
such that u G C'““'^(M”), hut u ^ C'^®’''"(i?i/ 2 )- Moreover, the same happens for the 
operators 01. 21) - 01. 4p . 

Thus, even when / = 0, solutions u to 

( Lu = 0 in i?i 
I u = g in ML \ Bi 

are in general no better than C^®(i?i/ 2 ) if g is not better than L°°(M”). This means 
that the interior regularity of solutions to 01.11) depends on the regularity of g in 
R”' \ when the kernels are not regular. 

Finally notice that, even in case that g is in all of R"" \ 12, this is not enough 
to deduce that u is regular inside 12! Indeed, because of 06.3p . one has to control 
the term ||m||c“(R") in order to have a estimate inside 12. For this, it is not 

enough to have ||fi'|lc“(R"\o) < C, but one also needs to control the regularity of u 
across (212. In other words, the boundary regularity of u is needed. We will come 
back to this in Section [HI 

®The estimate in [55] is a much stronger result, which holds for fully nonlinear equations and 
also for x-dependent equations. 
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7. Boundary regularity 


In this section we stndy the bonndary regularity of solutions to 

f Lu = / in n 
1 M = 0 in 


(7.1) 


We will hrst look at the optimal Holder regularity of solutions u near dVL, to then 
see more hne results for the class of kernels (11.51) . 


Remark 7.1. Once fl7.ll) is well understood, the boundary regularity of solutions to 
(11.11) follows from the results for (17.11) . at least when the exterior data g in (II.Ij) is 
regular enough. 

Indeed, assume that m is a solution to fll.ip and that g G \ O). We may 

extend 5 ^ to a function in M” satisfying g G and consider u = u — g. Then 

u satishes h = 0 in M"' \ O, and Lu = f — Lg =: f in Q. Hence, u solves fl7.ip with 
u and / replaced by u and /, respectively. 


7.1. Holder regularity for u. As we saw in Section [5], given any stable operator 
fll.2D - fll.5p . the function 

Uo{x) = (1 - 

is an explicit solution to fl7.1l) in hi = Hi; see Lemma [5.41 
The function uq belongs to C'^(Hi), but 

Mo ^ for any e > 0 . 

This means that, even in the simplest case K{y) = one can not expect 

solutions to be better than 

For the class of kernels fll.5l) . using the explicit solution uq and similar barriers, 
it is possible to show that, when hi is solutions u satisfy 

\u\<Cd^ in n, (i(x) = dist(x, M" \ hi). 


Combining this bound with the interior estimates, one gets the following. 

Proposition 7.2 (Optimal Holder regularity, [IZ]). Let L be any operator of the 
form fll.2p - fll.5p . and H any bounded domain. Let f G L°°{Q), and u be the 
weak solution of fl7.1l) . Then, 

||ii|lc'>(o) ^ C'll/IU°°(u) 

for some constant C that depends only on n, s, Vt, A, and A. 


Furthermore, using similar barriers one can show: 


Lemma 7.3 (Hopf’s Lemma). Let L be any operator of the form fll.2p - fll.5l) . and O 
any bounded domain. Let u be any weak solution to fl7.1l) . with / > 0. Then, 
either 

u > cd^ in Ll for some c > 0 


or u = 0 in Vt. 
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Thus, the boundary regularity of solutions depends essentially on the construction 
of suitable barriers, and the behavior of these barriers near dQ depends on the class 
of kernels under consideration. 

When the kernels K{y) are of the form fll.Sp . we have seen that solutions behave 
like near the boundary dQ. However, for the class of kernels fll.dp . the only barriers 
that one can construct behave like near the boundary, for some 0 < cio < s. 
This means that for this class of operators one can only prove u G 


7.2. Regularity of n/d®. We have seen that, when 


K{y) 


aiy/\y\) 


(7.2) 


then all solutions u behave like near the boundary. For this class of kernels, much 
more can be said about the regularity of solutions near dQ. 

Indeed, the quotient u/d® is not only bounded, but it is also Holder continuous 
up to the boundary. 

This regularity of u/d^ yields in particular the existence of the limit 


u u{x) 

— (z) := hm , , , 
d^ d^{x) 

for all z G dVL. This function u/d^ on dVL plays sometimes the role that the nor¬ 
mal derivative du/dv plays in second order equations. For example, it appears in 
overdetermined problems m, integration by parts formulas | 1 H], and free boundary 
problems [T2] . 

The hrst proof of this result was given in [15] for the case K{y) = and 

more recently the result has been improved by Grubb in |29l |30] and by the author 
and Serra in [l 6 l 07]. These results may be summarized as follows. 


Theorem 7.4 (|17l 06] |29l |30]). Let s G (0,1), L be any operator of the form fll.2p - 
fll.Sp . H be any bounded domain, and u be any solution to fl7.ip . with f G 

Then, depending on the regularity of the function a in fl7.2p . we have the following. 

(i) If a is any measure and H is , then 

\\u/d^\\cs-e(Jl^ < C'||/||^oo(Q). 

(ii) If a E and Ll is for some small 7 > 0, then 

\\u/d'^\\cs+j(Q'^ < 

whenever s + j is not an integer and f G G'’'(H). 

(hi) If a E G is C°°, and f E then 

\\u/d^\\Qs+j(ji) < G||/|| C'r{Q) 

for all 7 G (0, 00) such that s + ■j is not an integer. 
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Part (i) corresponds to Theorem 1.2 in HU, part (ii) was established in | 16 ] in the 
more general context of fully nonlinear equations, and part (iii) was established in 
[211 El] for all pseudodifferential operators satisfying the /i-transmission property. 
(In the excepted cases of (iii), more information is given in [30] in terms of Holder- 
Zygmund spaces C^, when s + 7 is an integer.) 

As in the interior regularity estimates, the regularity of the kernel K affects the 
regularity of the solution u. In this direction, even if one assumes that is C°° and 
/ G C'°°(f2), the result in part (i) can not be improved if the kernels are singular; 
see [i7] . 

On the other hand, notice that these results yield {u — g)/d^ G C'®’'''^(f2) for 
solutions u to fll.ip : see Remark l7.ll This gives a description of u near dfl up to 
order 2s + 7 . For example, if s + 7 < 1, this means that there exists a function 
b G O'^(cAl) such that 

|u(a;) — g{x) — b{z)d^{x)\ < C\x — for 2 ; G (90, x G Bi{z). 

When s + 7 G (1,2) one has a similar expansion with an additional term of order 
and more terms appear for higher values of 7 . 


8 . Further interior regularity 


As we saw in the previous Section, solutions u to flT.ip are or (7“° up to the 
boundary when K{y) is of the form fll.4p or fll.Sp . Moreover, the same happens for 
solutions to fll.ip whenever the exterior data g is regular enough. 

Hence, this means that solutions u will be globally Holder continuous, and there¬ 
fore we may apply the estimate 

||m||c“+27Si/ 2) - ^(ll/llc“(o) + ll'w||c“(R")) (8.1) 

to get that solutions u are ( 7 ^^+" inside for some a > 0 —more precisely, a = s in 
case fll.Sp . and a = ao in case fll.4p . This is enough to ensure that, when / G (7", 
any weak solution is a classical solution, in the sense that the operator L can be 
evaluated pointwise. 


The natural question then is: are solutions more regular than this? 


For the class of kernels fll.Sp the interior estimate fIS.ip does not give more than 
u G C[Q^^(fl), and we know that (18.1 p is sharp. Still, if one considers the solution u 
to 


( 8 , 2 ) 


Lu = 1 in 
u = 0 in 

in a C°° domain one a priori does not know if solutions are more regular than 
This question has been answered by the author and Valdinoci in 


Theorem 8.1 ([19|). Let L be any operator of the form fll.2p - fll.5p . and 12 be any 
bounded domain. Then, 

(i) If a is any measure and 12 is convex, then u G C';o^^*~^(12) for all e > 0. 
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(ii) If a E and ^2 is then u E for all e > 0. 

(iii) There is a (nonconvex) C°° domain f2, and and an operator L of the form 

fll.2p - fll.5p . for which the solution u to (18■2p does not belong to for 

any e > 0. 

Thus, solutions have different regularity properties depending on the shape of f2. 
Indeed, for operators L with singular kernels, the convexity of the domain gives one 
more order of differentiability on the solution. Also, notice that when a G then 
we gain interior regularity on the solution by assuming that the domain is C^’^. 

For kernels K in the class fll.dp . one would expect solutions u to be C'l+ss+ao-e jj^ 
the interior of domains. 
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